Dynamical Quantum Error Correction of Unitary Operations with Bounded Controls 
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Dynamically corrected gates were recently introduced [Khodjasteh and Viola, Phys. Rev. Lett. 
102, 080501 (2009)] as a tool to achieve decoherence-protected quantum gates based on open-loop 
Hamiltonian engineering. Here, we lurther expand the framework of dynamical quantum error cor- 
rection, with emphasis on elucidating under what conditions decoherence suppression can be ensured 
while performing a generic target quantum gate, using only available bounded-strength control re- 
sources. Explicit constructions for physically relevant error models are detailed, including arbitrary 
linear decoherence and pure dephasing on qubits. The effectiveness of dynamically corrected gates in 
an illustrative non-Markovian spin-bath setting is investigated numerically, confirming the expected 
fidelity performance in a wide parameter range. Robutness against a class of systematic control 
errors is automatically incorporated in the perturbative error regime. 
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I. INTRODUCTION 

The technological advances promised by quantum in- 
formation processing (QIP) science require an unprece- 
dented level of control in characterizing, manipulating, 
and measuring quantum systems [1]. In reality, no quan- 
tum system may be perfectly isolated from its surround- 
ing environment and available means for control and mea- 
surement arc unavoidably limited. As a results, devia- 
tions from the intended dynamical evolution arise, which 
are collectively referred to as 'errors'. Active quantum er- 
ror correction (QEC) methods aim to detect and counter- 
act the effects of errors during both quantum storage and 
processing of information via non-trivial quantum gates 
[2]. Given limited control resources, accurate quantum 
computation (QC) on an arbitrarily large (scalable) QI 
processor is only possible if QEC is implemented fault- 
tolerantly, that is, more errors is removed than is gener- 
ated when operations are themselves imperfect. The the- 
ory of fault-tolerant QEC (FTQEC) guarantees that this 
is possible in principle, provided that the error per gate 
(EPG) is below a minimum threshold [3, 4, 5]. Unfortu- 
nately, while architectures which can tolerate EPGs com- 
parable with experimentally achieved values have been 
proposed [G, 7], the corresponding resource overheads re- 
main prohibitive for existing QIP test-beds. Thus, QEC 
strategies which can effectively reduce errors at the phys- 
ical level are both crucial for boosting control fidelities 
in current quantum devices and, in the long run, for en- 
abling fault-tolerant QC. 

While traditional QEC methods achieve information 
protection by resorting to suitable redundant encodings 
and non-unitary quantum operations, dynamical quan- 
tum error correction (DQEC) aims to counteract the ef- 
fects of decoherence in a quantum system in a purely 
'open-loop' fashion, that is, solely based on suitable 
unitary (Hamiltonian) control. Historically, dynamical 
procedures for effectively removing unwanted couplings 
and/or evolution by time-dependent 'coherent averaging' 



have decades of tradition in nuclear magnetic resonance 
(NMR) spectroscopy [<^], resulting in an unmatched de- 
gree of coherent control in liquid-state nuclear-spin qubits 
[9, 10]. In particular, dynamical decoupling (DD) se- 
quences in NMR are early examples of DQEC for a closed 
quantum system, specifically designed to 'time-suspend' 
the underlying spin evolution by achieving a net trivial 
(identity) propagator, or noop gate. In the QIP con- 
text, DD methods have evolved into a powerful and ver- 
satile framework for achieving decoherence suppression 
in non-Markovian open quantum systems in a variety of 
control settings [11, 12]. While recent theoretical devel- 
opments include the exploitation of randomized design 
[13, 14], and the construction of concatenated [15] and 
optimal [1 (i] DD sequences, the growing practical signifi- 
cance of DQEC methods is witnessed by the intense par- 
allel effort in the experimental QIP community. Follow- 
ing simple proof-of-principlc demonstrations in a single- 
photon polarization interferometer [17], DD in its sim- 
plest (so-called 'bang-bang') form has been successfully 
implemented by now in solid-state nuclear quadrupole 
qubits [j -], fuUerenc qubits [10], coupled electron- nuclear 
systems [20], flying polarization qubits [!'], and, most re- 
cently, trapped ions [22]. 

Although the advances mentioned above demonstrate 
the usefulness of DQEC toward removing unwanted de- 
coherence in between control pulses and achieving a ro- 
bust NOOP implementation, full-fledged DQEC demands 
that decoherence be removed both in between and dur- 
ing control pulses, while implementing a robust non- 
trivial (non-identity) quantum gate. Dynamically cor- 
rected gates (DCGs) were introduced in Ref. 23 pre- 
cisely to meet this goal. As opposed to readily avail- 
able 'primitive gates,' which effect a unitary operation 
on the system with a 'bare' EPG, each DCG consists of 
a suitably designated sequence of control inputs in such 
a way that the desired target transformation is effected 
with a net reduced EPG. DCGs are inspired by both 
composite-pulse techniques for reducing classical control 
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errors (in particular, so-called 'fully compensating' or 
'universal rotation' sequences) [24, 25, 26, 27, 28] as well 
as strongly-modulating pulses [29, 30, 31] for high-fidelity 
'soft' pulses in NMR, while differing from either strat- 
egy in important ways. Specifically, DGCs compensate 
for error evolution due to a dynamical quantum environ- 
ment over which neither quantitative knowledqe nor di- 
rect control is assumed. Although DCG constructions 
arc useful only in the short gating time/weak coupling 
limit and cannot counter stochastic control errors, a key 
advantage is that no extra qubits or measurement capa- 
bilities are needed, thanks to the underlying open-loop 
control design. On the one hand, DCGs relate to ongo- 
ing effort on improving pulse-shaping methods in order to 
better accommodate physical constraints and/or increase 
robustness against different errors. In particular, DD- 
inspired constructions for pulses capable to 'sclf-rcfocus' 
unwanted couplings have been recently investigated in 
[32, 33, 34]. These efforts nonetheless consider a specific 
control task (noop implementation) and/or do not incor- 
porate the effect of generic quantum errors. DCGs, on 
the other hand, are applicable to any desired unitary op- 
eration, and have the flexibility to handle a generic class 
of non-Markovian error models on qubits. Furthermore, 
DCGs are designed having high portability in mind, so 
that they can replace the already available gates when- 
ever appropriate (mild) compatibility conditions are met. 
Substituting primitive gates with DCGs can then signif- 
icantly lower physical EPGs, ultimately paving the way 
for more practical FTQEC architectures. 

In this paper, we expand our analysis of DCG con- 
structions with the twofold goal of describing in detail the 
mathematical groundwork underlying the results of [23] , 
as well as further demonstrating the flexibility and po- 
tential of DCGs for high-fldclity quantum control in open 
quantum systems. The content is organized as follows. 
In Sec. II, we formulate the DCG synthesis problem in a 
control-theoretic setting, and describe the mathematical 
assumptions and concepts that are relevant to the sub- 
sequent steps. In Sec. Ill, we review the fundamentals 
of DD methods based on realistic bounded-strength con- 
trols, within so-called Eulerian DD [^io]. In particular, wc 
re-interpret the latter as a special DCG aimed at quan- 
tum state preservation, and prove a 'No-Go' theorem that 
restricts the portability of 'black-box' DCGs beyond the 
NOOP gate. Sec. IV shows how to evade the No-Go the- 
orem and describes general DCG constructions in detail. 
Special attention is paid to elucidating the compatibility 
requirements that emerge for physically relevant qubit 
error models, including arbitrary linear dccohcrcncc and 
pure dephasing. In Sec. V, a concrete spin-bath deco- 
herence scenario is used as a case study to numerically 
validate DCGs and to gain insight into open-system and 
control features which affect their performance. In Sec. 
VI, we elaborate on the effects of additional error evo- 
lution induced by control non- idealities and/or internal 
unitary drifts. A summary of the main findings and open 
questions conclude in Sec. VII. 



II. CONTROL-THEORETIC SETTING FOR 
DYNAMICAL ERROR CORRECTION 

Consider an isolated quantum system S with an in- 
ternal Hamiltonian Hs which can be controlled through 
the application of a time-dependent Hamiltonian iJctri(i)- 
This control enables us to perform certain quantum tasks 
on S. For example, by turning off Hs completely, an arbi- 
trary state of S may be preserved. Complex control tasks 
such as unitary quantum gates and algorithms composed 
of such gates may also be performed by modulating Hg 
through Hctr\{t)- In reality, however, 5* interacts with an 
environment system (referred to as the 'bath' B hence- 
forth) via an interaction Hamiltonian Hsb- In the pres- 
ence of Hsb, a designated control sequence on S is not 
sufficient to accurately implement a desired evolution, 
due to non-unitary decoherence errors that are induced 
once the bath degrees of freedom arc traced out. 

Let, as usual, the initial joint state of S and B be sep- 
arable, with p'" = P™®Pb^ ^nd = |'0in)(^in|- Ideally, 
the action of a target unitary gate Q effected over a gat- 
ing interval [0, T] should result in a final state of the form 
pI'^(T) = QVsQ ® PB^iT), for some (irrelevant) final 
state p^^{T) of the bath. In contrast, the actual action 
of the gate in the presence of Hsb results in the combined 
joint state PssCT), and the system ends up in an erro- 
neous reduced state p|(T) = Ttb{ps"{T)) ^ Q^'p'sQ- 
While different measures may be envisioned to quantify 
the deviation between the intended and the actual evolu- 
tion, the resulting EPG will be proportional to the gating 
time T, as well as to the overall 'error strength,' appro- 
priately defined. In essence, DCG constructions aim to 
improve the scaling of physical EPGs with the gating time 
in the limit where the latter is sufficiently small, so that 
a perturbative approach is meaningful. That is, a gate 
Qdcg is an i-th order DCG if it obeys 
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(1) 



epgM 

where the asymptotic 0-notation is used to signify the 
limit where the zero-th order bare EPG is close to 0. 
If desired, the above defining condition for a DCG may 
be expressed in terms of gate fidelity improvement once 
a specific fidelity metric is chosen and is related to the 
EPG amplitude in Eq. (1). Let, for instance, fg denote 
the worst-case fidelity in implementing Q [2] , 



fo = min Tr- 



(2) 



By taking the 'infidelity' (1 — /g) as a measure of the 
gate error probability, a ^-th order DCG implies an im- 
provement ratio on the order of [13, 2G, 3G] 



Wdcg 



(3) 



where /qdcq is the fidelity of implementation of a se- 
quence of primitive gates whose fidelities are all similar 
and given by /q. 
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In this paper wc will only focus on first-order DCGs 
{£ = 1), by postponing the construction and character- 
ization of higher-order DCGs to a separate forthcoming 
analysis [37]. 

A. Error model assumptions 

Let iS* consist of n qubits, and let B represent another 
quantum system over which knowledge is limited and no 
direct control is available. The free evolution of S and B 
in the combined Hilbert space H = Hs^Hb is described 
by a bare internal Hamiltonian Hint of the form 

Hint = Hs®lB + HsB+Is®HB, (4) 

where Is (Ib) is the identity operator in the operator 
space B{Hs) [B{Hb)], respectively. Otherwise obvious 
tensor products components will be dropped from now 
on; e.g., Hb will be understood as Is ® Hb- The inter- 
action Hamiltonian Hsb is responsible for decoherence 
effects whereby initially separable states of 5* and B be- 
come entangled and loose their capacity for QIP. The 
typical time scale over which such effects become signif- 
icant (loosely referred to as the 'decoherence time') de- 
pends in general on both Hs b and various features of the 
bath, including p'g [38, 39]. Without loss of generality, 
HsB can be expanded as Hsb = J2-f ^-y ® -^7j where 
{Sj} is an (Hermitian) operator basis for BiJis), such as 
arbitrary products of Pauli operators on the i-th qubit, 
denoted by {S'i^} (a e {I,X,Y,Z} = {0,1,2,3}), and 
B^ being non-zero but otherwise unspecified bath opera- 
tors. In most circumstances, the dominant terms in Hsb 
belong to a restricted subspace of all possible coupling 
operators. In particular, the linear decoherence model is 
defined by letting: 

n 3 

^iB = EE^"^®^"^- (5) 

2—1 O— 1 

Physically, the class of error models described by Eq. 
(5) allows for arbitrary combinations of phase- and 
amplitude- damping processes involving a generic axis 
and degree of correlations between different qubits (from 
independent to fully collective linear interactions) [40]. 
In DQEC approaches, it is useful to associate an error 
subspace fig C BiTL) to the set of error generators wc 
wish to suppress. For linear decoherence, 

fie =span{5W ,^B|Va,i; B G 6(Hb)} = ^^^'^ (6) 

includes arbitrary system-bath operators that induce 
single-qubit errors on S. Also notice that the strength of 
the system-bath interaction Hamiltonian is constrained 
to be sufficiently weak in order for both QEC and DQEC 
to be useful in practice. This motivates assuming that 
both Hsb and Hb arc bounded in an appropriate norm, 
although otherwise potentially unknown. The unitarily 



invariant operator norm [ ' "I] will be used throughout 
our analysis. 

If Hsb is not present, a desired unitary operation Q on 
S can be generated, as mentioned earlier, by adjoining a 
semi-classical controller, described by a time-dependent 
Hamiltonian Hcu\{t). Depending on the system and con- 
trol specifications, certain components of the internal 
system Hamiltonian might be essential to ensure univer- 
sal controllability, whereas other will induce unintended 
(unitary) error evolution. In general, we may thus de- 
compose Hs = Hs.g + Hs^ei where the separation into 
the gating and error contribution depends on the target 
transformation Q; for instance, Hs = Hs.e for noop, 
Q = Is- In the presence of Hsb, let the total error 
Hamiltonian be defined by 

He = Hs.e + Hsb + Hb, (7) 

and, correspondingly, we represent the control action in 
terms of a gating Hamiltonian of the form 

TJgatc (t) = iJgatc (t) (^Ib^ i^ctrl (t) + Hs,g. (8) 

In principle, a perfect implementation of Q could still 
be ensured through infinitely strong, instantaneous gat- 
ing Hamiltonians, during which H^ has no time to act. 
However, this ideal scenario cannot be achieved in real- 
ity, although it may be realized approximately. Even if 
Hsb 7^ 0, it may be useful to have a way for compar- 
ing operators in B{'H) which differ only through directly 
uncontrollable pure-bath terms such as Hb. This can be 
accommodated by defining equivalence classes of opera- 
tors: We say that X is equal to Y modulo the bath (modi? 
for short) iff their difference is a pure-bath operator, 

X = Y mods 3Xb : X-Y = Is®Xb, 

where Xb is an operator acting on B. For example, 
Hgateit) + Hb = Hgate{t) modi?. With this in mind, 
we use the words ideal or desired to refer to cases with 
He = mods, where gates can be perfectly implemented 
by a suitable closed-system protocol Hgi^toit) over time 
T. Furthermore, we will primarily focus on open-system 
decoherence errors, that is, we will treat Hsb as the lead- 
ing source of errors in Eq. (7), by assuming that: 

(al) Perfect control assumption: No additional 
errors are introduced by the controller; 

(a2) Driftless system assumption: No additional 
errors are introduced by the system's internal evolution. 

Assumption (al) implies that the applied control 
Hamiltonian i?ctri(i) is perfect, subject to operational 
constraints to be specified later (Sec. II B). While re- 
alistic controls will suffer in general of both systematic 
and stochastic (random) imperfections, these errors are 
intrinsically classical in nature. Likewise, assumption 
(a2) is partly justified by the fact that Hs is fully known 
and, in the closed-system limit, complete control over S is 
available. Once DCG constructions arc found under the 
above simplifying assumptions, it is at least in principle 
conceivable that simultaneous compensation of classical 
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and quantum (unitary and decoherence) errors could be 
achievable by suitably merging DCG constructions with 
existing robust control techniques, although several non- 
trivial aspects may remain. While additional discussion 
is deferred to Sec. VI), the impact of control imperfec- 
tions will be also numerically assessed in Sec. V. 



B. Control assumptions and error measure 

Our goal is to use available control resources to gener- 
ate unitary gates on S while minimizing the sensitivity of 
the gate action to the error Hamiltonian H,,. In practice, 
allowed manipulations will be restricted by a number of 
constraints, which can be identified by writing 



(9) 



in terms of a subset of admissible control Hamiltonians 
{Vu} and of control inputs {hu{t)}, the latter being real- 
valued controllable functions of time. Limited 'pulse- 
shaping' capabilities may additionally restrict the val- 
ues attainable by the control inputs hu{t) and/or their 
derivatives. Two main constraints will be incorporated 
into the present DCG design: 

(cl) Finite power, which requires bounded control 
amplitudes, hu{t) < /imax; 

(c2) Finite bandwidth, which requires a bounded 
Fourier spectrum, hence a minimum time scale Tmin for 
modulation. 

A 'primitive gate' will refer to a physically available 
gate used to generate more complex transformations. In 
the simplest case, a primitive gate may be implemented 
by turning on and off a single control Hamiltonian in 
the available set, subject to the above assumptions. A 
'control block' will likewise refer to a time interval dur- 
ing which multiple primitive gates are applied sequen- 
tially in a predetermined manner. From the standpoint 
of universal QC, two different scenarios arise depending 
on whether the set of switchable control Hamiltonians 
{Vu} in Eq. (9) supports a universal set of single- and 
two- quhit gates. If this is the case, {Kt} generates a 
dense subset in SU(2"), thus arbitrary unitary gates on 
S may (ideally) be approximated to a desired accuracy 
by turning on and off primitive Hamiltonians in the set. 
Note that if no contribution from the internal Hamilto- 
nian Hs is needed to achieve universality {Hs,g = 0), 
then to the purpose of DCG design we may effectively 
treat Hs as a pure error (Hs = Hs.e) and the problem 
as effectively driftless, in line with (a2). 

Let Q denote a quantum gate generated through the 
application of Hcti\{t) over the time interval [0, T]. Re- 
calling Eq. (8), the gating propagator J7gato(i,0), t £ 
[0,T], can be defined in general as 



C/gate(i,0) = T+ 



cxp 



Hgi,teis)ds 



such that f/gatc(7', 0) = Q. We use Ugate{t) to denote 
Legate (t, 0) when there is no ambiguity about the starting 
time of the gate. The overall propagator for the system 
plus bath is given by 



c/(r,o) = T+ 



cxp ( - I I (Hgatois) -\- He)ds 



A natural way to isolate the effect of the error Hamilto- 
nian is to effect a transformation to a frame that 'toggles' 
with the applied control, that is, to move to the interac- 
tion picture generated by J7gatc(i, 0): 



C/(i,0) 

-i$(t,0) 



C/gate(t,0)e-^*(*^°), 



cxp 



I I Heis,0)ds 

'0 



(10) 
(11) 



where the 'modulated' error Hamiltonian reads 



i/e(t,0) = C/gatc(t,0)ti7eC/gatc(t,0). (12) 

The Hermitian 'error action operator' ^(T, 0) defined in 
Eq. (11) is directly related to the effective Hamiltonian 
which describes the joint evolution in the toggling frame. 
Mathematically, ^{T, 0) provides, up to pure bath terms, 
a measure of the overall deviation of the actual evolution 
of S from the desired unitary evolution Q. Since the 
norm of ^{T, 0) may be used to lower-bound the mini- 
mum achievable fidelity [42] , the latter may be used as a 
natural EPG metric for DCG constructions [Eq. (1)], and 
the associated quantum-control problem may be viewed 
as the minimization of $(r, 0), up to purc-bath terms. 
The solution obtained in [2:!] and further investigated 
here is perturbative and analytic: We guarantee that 
$(7,0) ^ EPG'^' is of second order in Tmin, and refer 
to this perturbative cancellation as 'correcting first-order 
errors,' while similarly referring to the residual (asymp- 
totically smaller) terms in $ as 'uncorrected'. 



C. Tools for error analysis 

Our next step is to describe how the error action <& 
can be evaluated and approximated for complex evolu- 
tions involving a composite control block. Suppose that 
the interval [0, T] is decomposed into two segments [0, Ti] 
and [Ti , T] . By using the composition property of uni- 
tary propagators and, for each interval, the separation 
between intended and error contribution, we have 



-i*(T,0) 



Ul,,{T, 0)C/gatc(r, T,)e-^^™U{T,,0) 

exp{-i C/g^ate,l*2t/gatc,l) exp(-i$i), (13) 



where 



gate,l 

-2^1 _ 



C/gate(Ti,0), 
„-j*(Ti,0) ^■ 



-i*(T,Ti) 
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Eq. (13) allows to recursively combine segments of evo- 
lution entering a composite unitary gate and to calculate 
the EPG associated with the combined sequence in terms 
of EPG of each segment. Physically, this is equivalent 
to transforming to consecutive toggling frames. In this 
sense, the error associated with each individual segment 
can be isolated and computed, provided that the con- 
trol propagator during the segment and the control path 
applied up to the beginning of the segment are known. 

Generally, when N gates Qi, • • • , Qn are applied in 
succession, the error operator associated with the combi- 
nation A = QnQn-1 ■ ■ ■ Qi is given by with 

exp(-i$A) = e-*'^"-i*««-f'"-i ••• (14) 

^-^P!'^>Q^Pl^-^P^l^^>Q,Po^ 

where the 'partial' gating propagators up to segment i 
along the path are given by 

P. = Q^Q^~l • • • Qi, i = 1, . . . , AT - 1; Pa = Is- 

Notice that from Eq. (14), in the limit of infinitesimal 
segments, we recover the toggling frame expression for 
$(r,0), as given in Eqs. (11)-(12). Eq. (14) also lends 
itself naturally to approximation. If the errors associ- 
ated with the individual gates Qi arc known, then the 
(discrete-time) Magnus expansion provides a systematic 
(albeit computationally demanding for higher orders) se- 
ries expansion for <I>^: 

oo 

i=l 
N 

$11' = Y,p.U<^Q^p^-i, (15) 

oo 

= 0(max|l$Q.f), 

i=2 

and the above estimate for higher-order corrections [43] 
holds as long as (absolute) convergence is ensured, that 
is, if A^max|j$Qj| < tt. The (continuous-time) Magnus 
expansion can similarly be applied to approximate errors 
directly in the toggling frame [12, 35] or to estimate the 
individual EPGs ^q^- Starting with Eq. (11), in analogy 
to Eqs. (15) we have: 

oo 

$(t,0) = 5] $[''(*, 0), 

$W(t,0) = / He{s,0)ds, (16) 
Jo 

\\Y.<f^Hm\\ = o(iiti/eii'), 

1=2 

as long as t \\He\\ < tt. In what follows, wc will use Xl"! to 
denote the n-th Magnus expansion term for the operator 
X and to denote the sum of n-th and higher-order 



terms. Furthermore, when depicting the control inputs 
for a gate, we will assume that the gating period starts at 
the time to = 0, unless otherwise stated. This is possible 
due to the fact that when individual gates are consid- 
ered, the EPG does not depend on the initial gating time 
unless the physical error Hamiltonian in Eq. (7) is itself 
explicitly time-dependent. 

III. DYNAMICALLY CORRECTED NOOP 

The 'no operation' (noop) gate plays an especially im- 
portant role from both a quantum-control and a QIP 
standpoint. On the one hand, even in a closed-system 
setting, the ability to un-do the effect of the drift and 
realize an arbitrary closed control trajectory is a pre- 
requisite for complete controllability [44]. On the other 
hand, preserving an arbitrary quantum state on all or 
part of a quantum circuit for a desired storage time is 
an essential subroutine for complex quantum algorithms. 
As such, achieving the NOOP gate is an important bench- 
mark for DD methods. In turn, DD provides us with 
an inner combinatorial structure for DCGs which goes 
beyond the NOOP gate itself. 

A. Basics of dynamical decoupling 

DD is the most prevalent DQEC strategy. In essence, 
DD schemes subject the system to either sequences of in- 
finitely strong instantaneous 'J-pulses' (bang-bang DD, 
BB DD) or bounded-strength control modulations (Eu- 
Icrian DD, EDD), in such a way that undcsired evolution 
induced by the bath is coherently averaged out over a 
time scale longer than a minimum control time scale Tc- 
While DD performance can depend heavily upon the de- 
tails of the sequence used, a basic requirement for averag- 
ing is that Tc be sufficiently short as compared to the time 
scales of the erroneous dynamics to be removed. While, 
as mentioned, a myriad of progressively more elaborated 
DD schemes have emerged [13, 45, 46, 47], we briefly re- 
view here the original group-theoretic construction [12] 
in a form suitable to our current scope. 

Consider a finite group Q = {gi}f^i, d = \Q\, which 
acts through a faithful unitary (projective) representa- 
tion {Gi}f^i on Tis and through a corresponding adjoint 
representation on B{Ti.s)- These representations can be 
extended trivially to the combined Hilbert space Ti by 
using {Gi (8) Ib} as the representation of G- 5 is a de- 
coupling group for a subspace Q C B{Hs) of (traceless) 
operators on Tig iff 

d 

VEen-. ng{E) = Y^ g\eg, = o, 

i=l 

where ^ng(-) is the projection superoperator onto the 
tj-invariant sector, and we refer to each Gi as a decou- 
pling operator for Vl. For example, consider = Z2, 
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represented by {I,X} on the state space Hs = of 
one qubit. This group is a decouphng group for the 
subspace nj_x = spa,n{Z,Y}, thus it can correct er- 
ror Hamiltonians in the correctable error space flc = 
{Z (E) Bz + Y ^ By + Is ^ Bq}, regardless of the arbitrary 
tensor components Ba- In general, if fig is the subspace 
of system-bath operators generated by all errors wc wish 
to remove, a good DD group ensures that 3 ^^e- 

Consider next a control propagator Ubd (t) applied on 
d consecutive time intervals li = [ti-i,ti\, where i = 
to = 0, and (t^ — ^i-i) = t, such that each 
decoupling operator is implemented once, beginning from 
the identity Is = Gi. That is, 

In principle, the discontinuous (in time) propagator 
C^DD(i) can be realized by a singular control Hamilto- 
nian that is switched on at ti with infinite strength, and 
instantaneously rotates S according to Di = Gi+iG\ for 
i = l,...,d— 1, and Dd — G\. Notice that the propaga- 
tor after td is the identity operator, which amounts to a 
NOOP gate associated with the whole evolution. S and B 
undergo free evolution governed by during {ti-i,ti) 
with an error given by but S is kicked at ti by Di. 
Let T = dr. The error associated with the combined 
sequence in Eq. (17) can be obtained by using Eq. (11): 

d 

= n exp(-zTGli/eG,0 = <i>W^ + 

i=l 

$W = exp(-zrX:GlM.)=e-^n-(«^),(18) 
1=1 

where the Magnus expansion of Eq. (15) has been used 
to approximate in the leading powers of tH^, as- 

suming that T||7Je|| < TT. As long as € Jig, then 
*^DD ~ = modi?, thus the protocol will produce 
a (first-order) error-corrected NOOP on S. Notice that 
we could alternatively envision the DD sequence as a se- 
quence of (ideal) gates Di with errors = (mod B) 
and primitive NOOP gates F (consisting of free evolution) 
with error = r_ffe, apphed as DdF ■ ■ ■ D2FD1F, with 
a combined error given by Eq. (14). The uncorrected 
error associated with the DD sequence is given by the 
higher-order Magnus terms [12, 48]. 

Theoretically, the cfhciency of the basic 'periodic' BB 
DD scheme based on (17) is limited by the requirements 
that (i) the control time scale Tc = dr be short enough 

for the Magnus expansion to be applicable; and (ii) the 

[2+1 

higher-order error terms 'I'dd be sufficiently small. The 
impact of higher-order corrections is especially impor- 
tant if repeated control cycles are implemented, in order 
to achieve a net NOOP over a desired (arbitrarily long) 
storage time. More sophisticated DD techniques address 



these issues by invoking a number of design improvements 
- including recursive constructions [45] , randomization of 
the control path [13, 49], and optimization of the control 
intervals [4G, 47]. Presently, however, for the generic lin- 
ear decoherence model of Eq. (5), no provably optimal 
DD sequence exists; that is, unlike in the pure dephas- 
ing case where decoherence occurs in a preferred basis, 
no efficient high-order DD schemes are known. Further- 
more, in complex systems of interacting qubits, DD may 
be used to remove unwanted inter-qubit couplings or to 
selectively enforce a desired dynamics. In such scenarios, 
identifying a minimal DD group and its corresponding 
physical representation is, likewise, a largely open prob- 
lem in combinatorics [50]. 



B. Eulerian decoupling design 

A major simplification afforded by the BB setting is the 
separation between the evolution due to the controller 
and the one due to the internal open-system Hamilto- 
nian i?int- Realistic operations, however, always entail 
a bounded control amplitude, thereby a finite duration. 
Even if devices exist where 5-pulses are an adequate first 
approximation, from a control-theory standpoint it is 
highly desirable to recover the BB setting as a limiting 
case of a formulation based on (non-singular) admissible 
controls. This is the central motivation in EDD [35]. 

Consider a set of generators F = {hjYjLi for Q, such 
that each element g (z G can be written as an ordered 
product of hi's. The Cayley graph G{G,T) of Q with 
respect to F pictorially represents the generation of all 
elements of Q through application of hj : Each vertex of 
G(fJ,F) represents a group element, and a vertex is 
connected to another vertex gg' by a directed edge 'col- 
ored' (labeled) with the generator h iff gi' = hgg. The 
number of edges of G(C/,F) is thus given by L = dm. It 
can be shown that every Cayley graph possesses an Eu- 
lerian cycle that visits every edge exactly once. Let C be 
an Eulerian cycle on G{Q, F) that starts (and ends) at 
the identity (see Fig. 1 for a pictorial view). 

This cycle can be associated to the control propaga- 
tor C/EDD(t) which describes evolution in L consecutive 
intervals Lj = [tj_i,tj], that is, 

%DD(tj) =^jt^EDD(tj-l), J = l,...,i, (19) 

where Fj is the group representation of the generator 
labeling the j-th edge in C. According to Eq. (19), the 
control propagator follows the Eulerian path faithfully: 
The gating Hamiltonian i?gatc(t) varies in such a way 
that each evolution segment labeled by a generator hj 
along the path is implemented by a corresponding gate 
Fj. The evolution over Tc = Lt thus consists of each 
generator Fj £ F being implemented precisely d times. 
Suppose that the error associated with each gate Fj is 
given by $ , determined by the manner in which Fj is 
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Figure 1: (color online) The Cayley graph G{Q,T) with 
g = Z2®Z2 = {9i}U = {(0,0), (0,1), (1,0), (1,1)} and 
r = {/ii,/i2} = {92,93} = {(0, 1), (1,0)} with an Eulerian 
path of length L = 8 depicted. This group with the Cay- 
ley graph and the Eulerian cycle depicted here will be used 
in Sec. IV C to implement DCG constructions for the linear 
decoherence model. 



physically implemented: 

exp(-i<I>F_,) = r+ CXp(^-iJ^ C/gato,j(0^-ffe[/gatcj(i) 

where we have assumed that Fj is generated through the 
application of a gating Hamiltonian i/gate,_F, (t) in the in- 
terval [0, r] (which can be shifted forward in time), with 



Ugatc, j{t) — T + 



exp 



'0 



The combined error associated with the EDD sequence 
can be approximated by a Magnus expansion [Eq. (15)]: 



(20) 



with 



EDD 



0(maxj ||$i.J ). In EDD, we require 

that the errors $ p. associated with the control generators 
belong to the error set fie corrected hy Q. As long as 
this is the case, then $edd ~ ^ modS, and the EDD 
sequence performs a NOOP with an asymptotically small 

[2-1-1 

error $edd compared to the uncorrected evolution. 

Given arbitrary linear decoherence, we will show in 
Sec. IV C how readily available control options for real- 
izing the generators of G can produce errors that indeed 
belong to Qe- For a generic open system, it need not 
be apparent which type of Hamiltonians and control in- 
puts (if any) might be used to generate EDD generators 
such that their respective errors arc correctable by 
(some) G, a question we refer to as the 'gate permissibil- 
ity problem'. To summarize, provided a set of permissible 
generator gates can be identified, EDD yields a construc- 
tive solution to our main problem as defined in Sec. II B, 
with the desired target gate being the NOOP gate. 



C. No-go theorem for black-box constructions 

An important feature of EDD is that the net error 
given by Eq. (20) is canceled irrespective of how each gate 
Fj representing a decoupling generator is implemented in 
terms of the control inputs, provided that (i) the same 
implementation is used each time hj appears along the 
path and (ii) each error is correctable by Q. To state 
it differently, no a priori relationship between the errors 
{^Fj} 'is assumed, as they cancel independently up to 
the first order [Eq. (20)]. In this sense, error correction 
in EDD is oblivious to the details of the control. Unfor- 
tunately, as we shall prove next, such 'control-oblivious 
decoupling' can only dynamically correct the NOOP gate. 

Consider, specifically, a control sequence composed of 
N gates Qi applied on S in the presence of He- Ideally, 
let the sequence be intended to implement a gate A = 
Qn ■ ■ ■ Qi- The EPG <I>q. associated with Qi is assumed 
to be an arbitrary function of the gate Qf. 



Qi 



£{Q^). 



£ depends in principle on as well as the details of im- 
plementation. Suppose, however, that as in EDD no in- 
formation about such dependence is directly or indirectly 
incorporated into control design. We have the following: 

Theorem. Consider a sequence of gates A = 
QnQn-i ■ ■ ■ Qi, o.nd suppose that the combined first- 
order error $^ for this sequence is equal to zero modB 
as long as the individual errors ^q- belong to a subspace 
rig C rig of dynamically correctable operators for all i. If 
no further assumptions are made on $q . , then the expec- 
tation values of operators in $7^ are preserved by A. 

Proof. Since no assumptions are made on (hence 
on (^Q. modi?), the function £ from unitary system oper- 
ators into fie is arbitrary. In particular, we may consider 
two error models defined by 

£i{U)^X, £2{U)^U^XU, 

where X is any operator in Vt^- Let ^a,i ('^'^,2) denote 
the combined error obtained under the error function £1 
{£2)- Let us define Pi = Qi - ■ ■ Qi for N > i > and 
Pq ~ Is- From Eq. (15) and using the fact that Qi — 
PiP}_i, we have: 

N 

= = (21) 

i=i 

N 

i=l 
N 

= P}XP, = mods. (22) 



Comparing Eq. (21) and Eq. (22) yields: 

VX e f7e : P^XPa - pIjXPn = AT - A^XA = modS. 
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We can write X a,s X = X + Is Bx where TrgX = 
so that X has no pure-bath component. The operator A 
acts on the system only, thus we have 

X-{A(x)lByx{A®lB)^OinodB =>X - A^^XA^O 

=>X - A'^XA = 

Clearly, A® Ig acts as identity on fig and can generate 
no non-trivial dynamics on this subspace. I 

In other words, we may yet effect a control-oblivious 
unitary transformation on the system but it will neces- 
sarily commute with the subspace of errors that Q can 
correct for, and will thus coincide with identity (noop) 
if the latter acts irreducibly on Tis- Nonetheless, one 
might be interested in suppressing errors from a partic- 
ular subset of more general errors. The No-Go theorem 
then still allows such a black-box solution to be viable 
for constructing primitive gates, as well as it allows to 
suitably combining DD with quantum encoding. For ex- 
amples of constructions where decoherence is suppressed 
using BB pulses and the desired computational dynamics 
is encoded in a subspace or subsystem see [43, 51, 52, 53]. 

From a practical standpoint, the No-Go theorem places 
restrictions on the ability to making fully portable DCGs 
where a fixed sequence of arbitrary implementations of 
gates effects a generic unitary gate on the system. In 
order to go beyond NOOP in constructing DCGs, we thus 
require to include certain specified gate implementations 
that circumvent the arbitrary EPG assumption. 



IV. DYNAMICALLY CORRECTED GATES 
BEYOND NOOP 

A central point in EDD is that the combined EPG is 
determined by the variation of the control propagator 
between the gates at the vertexes of the Cayley graph 
underlying the scheme. In view of the above No-Go The- 
orem, identifying a procedure where an analogous error 
cancellation is achieved but a non-trivial target unitary 
is effected on S*, requires access to relationships between 
errors associated to different gates. Identifying and ex- 
ploiting such relationship is possible because both the 
algebraic structure of i?e and the gating Hamiltonian are 
known. Thus, while the No-Go does not pose a funda- 
mental obstacle for DQEC, it does complicate the design 
of robust control schemes and reduces their portability. 
It is the goal of this Section to provide explicit proce- 
dures for constructing and analyzing DCGs, by first il- 
lustrating the general approach and obtaining estimates 
of the uncorrected (higher than first-order) net error, 
and next providing full detail on linear decoherence- and 
dephasing- protected DCGs. Both for added clarity, and 
because the presence of a non-trivial Hs 7^ requires a 
system-dependent analysis, we take advantage in what 
follows of assumption (a2), that is, we explicitly let 

Hs = 0, He — HsB + Hb, i?gatc(i) = -ffctri(i), (23) 



and further discuss the role of Hs in DQEC constructions 
in Sec. VI. 



A. DCG constructions and error estimates 

Probably the simplest (albeit not unique) way to en- 
force the constraint of non-trivial gate error relationships 
is to imagine that two different gates share the same EPG 
to lowest order. Let Q be the intended (primitive) tar- 
get gate, with associated error ^q, and imagine that a 
special NOOP gate Iq is available, such that the error 
associated with Iq is also $q: 

Uq^Q exp(-i$Q), Iq = Is cxp(-i$Q). (24) 

That the above requirement can in fact be fulfilled will be 
addressed in Sec. IV B. Assuming for now that Eq. (24) 
holds, consider attaching a self-directed edge Iq to each 
vertex on the Cayley graph G(C/, F), and let the Eulerian 
path C be modified so that the new edges added at each 
vertex are incorporated. The error associated with this 
new sequence (which clearly still acts as noop) is given, 
up to the first order, by 

*l^LG = *EL + ne($Q), 

and is canceled (modS) as long as the errors <I>q and ^p. 
belong to fig. Next, we change the final gate/edge Iq in 
the path we just created with the gate/edge Q with error 
<f>Q. By construction, the error associated with the com- 
bined sequence remains unchanged, but the sequence now 
implements Q as opposed to Is- Thus, we have succeeded 
in obtaining a DCG that performs the desired gate Q on 
S with a smaller error as compared to its uncorrected 
implementation. 

In line with the philosophy underlying DQEC ap- 
proaches, the errors due to Hf, are suppressed asymp- 
totically. This requires that the EPG of the primitive 
(uncorrected) gates be small to begin with in order to ex- 
pect a relative improvement in DCG performance [recall 
Eq. (1)]. Even if this condition is met, the uncorrected 
errors associated with a large number of operations per- 
formed in succession tend to inevitably build up, thus ad- 
ditional precautions must be taken to avoid catastrophic 
error growth. Within FTQEC, this can still be prevented 
provided that DCG residual errors remain correctable by 
the given code architecture and smaller than the relevant 
threshold. This prompts us to estimate the uncorrected 
errors associated with DQEC constructions for each gate. 

The asymptotic nature of DCG constructions rests on 
the validity of the Magnus expansion as an approxima- 
tion of the effective Hamiltonian for short time intervals. 
When used to estimate the EPG associated with a gate 
Q of duration T = Nt, a sufficient condition is 

\\H,^\T<n, (25) 

where A'' = 1 corresponds to a single primitive gate seg- 
ment. Two remarks are in order. First, in view of as- 
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sumption (c2), Tmin effectively limits the strength of er- 
ror Hamiltonians that can be corrected using DCGs. Sec- 
ond, while including the environment Hamiltonian Hb in 
He has the advantage of simplicity, the resulting conver- 
gence radius (time) for the Magnus expansion tends to be 
overly pessimistic in this way, in particular it may grossly 
overestimate the error in the case where \\Hb\\ ^ ||i/5B||. 
One solution is to move to a toggling frame that removes 
the evolution of B. We refer for instance to [4(S] for an 
analysis along this lines, yielding a tighter convergence 
radius for the Magnus expansion, Iji/gsll T < tt. 

The uncorrected errors that do not cancel as a result 
of a DQEC protocol can be bounded by the second- and 
higher- order contributions to the Magnus series [54] . Re- 
calling Eq. (12), we can estimate the uncorrected error 
^Q.DCG associated with a DCG Q as 



Q,DCG 



1 
2 



T pt 



[He{ti),He{t2)]dhdt2 



< — max||[i/,(ti),He(ti)]|| (26) 

4 ti<r2 



< 



2 WHbW \\H: 



1SB\ 



\H. 



SB I 



where we have used the unitary invariance properties of 
||-||. While Eq. (26) gives an upper bound for the uncor- 
rected error associated with any gate Q for which = 
modi?, it fails to capture many interesting features of the 
interplay between the pure environment evolution and 
the system-bath coupling. We will reconsider the latter 
numerically in the specific setting of Sec. V. 

Note that the first order error <i>g associated with a 
primitive gate Q (not corrected through DQEC) will nat- 
urally depend on the number of qubits n in the system. 
For the linear decoherence model we have 



1> 



(X 71, 



thus the EPG per qubit is expected to be constant (size- 
independent) up to the first order. For DCGs, the first 
order error is zero, and to obtain a dependence on the 

[21 

n, we need to focus on dcq, estimated above. Using 
Eq. (26), we may write 



[2+] 



,DCG 



Oin"), l<s<2. 



The upper bound is given by ||i?ss|j [Eq. (26)], whereas 
the lower bound is achieved, for example, when the bath 
operators Ba'' act locally, [Ba\ B^^''] cx Sij, for all a,f3. 



B. Finding gates with same error 

As described in Sec. IV A, a fundamental ingredient 
in DCGs is that two control sequences are found for ev- 
ery desired gate Q: One designated to perform a noop 



and the other designated to perform Q, such that the 
corresponding errors are equal up to the first order. We 
now give a simple method for generating such sequences, 
which in the language of Ref. '^7 realize a so-called first- 
order balance pair. 

Suppose that a given control input hf^it) is intended 
to produce a unitary gate Q during an interval [0,r], 
and let iJ^jp, and C/gatc(0 denote the gating Hamilto- 
nian and corresponding propagator for Q over t G [0, r], 
respectively. The error associated with Q is given by 



exp(-i$Q) = 



cxp I — I 



He{s,0)ds 



with He{t, 0) being the toggling frame error Hamiltonian 
as in Eq. (12). We next smoothly extend the control 
profile for Q to [0, 2t] by applying a gate Q' immediately 
after Q, in such a way that the overall control propagator 
for t € [0, 2r] is given by: 



_ / C^gatc(i), < i < r; gate Q, 

C/gatc(2r - <), r < < < 2t; gate Q' . 



The extended (composite) gate Iq = Q'Q implements 
the identity (noop), with a first-order error given by: 



*Q'Q - / C/U(t,0)tHeC/^atc(i,0)dt (27) 

-2/ C/gate(0^^eC/gatc(Odi = 2$W 
JO 

Consider next another control propagator associated 
with a gate Q1/2 over the interval [0, 2t] obtained by 
letting 

(^gate(i,0) = [/gate(t/2), < t < 2r. 

The first-order error associated with Q1/2 reads 



2t 



[1] 

Ql/2 



C/;,t,(i,0)tHeC/^'atc(i,0)dt (28) 



= 2/ C/gatc(i)tiieC/gatc(Odi = 2$gl. 

Jo 



By comparing Eqs. (27)-(28), we observe that while Q'Q 
ideally implements the identity and Q1/2 implements Q, 
the EPGs resulting from ^ arc equal up to the first 
order. 

Given enough control over the {ha(t)}, simple recipes 
may be given for implementing both Q' and Qi/2- For 
instance, Q' may be reahzed by applying —h^{t — T) (the 
reverse anti-symmetric profile to h'^{t)) when t > t, 
see Fig. 2. Likewise, to implement Q\j2i we may 

use a control input (t) such that /iq (s)c?s = 

J^^"^ h'2{s)ds. For example, if the pulse profiles are rect- 
angular, Qi/2 may be generated by using a control in- 
put with half the speed and half the strength of the one 
used for Q. The two prescriptions so obtained for QQ' 
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Figure 2: (color online) Control profiles for implementing two 
gates Iq — Q'Q and Q1/2, that correspond to a target identity 
gate and a target gate Q, with equal first-order error. 



and Qi/2 corresponds then, respectively, to the piecewise- 
constant profiles given in Eqs. (4)-(5) of [2->]. If we as- 
sume that all primitive gates have the same duration r, 
then the number of primitive gates required in each DCG 
will be dm + 2d^ d{m + 2) [recall Sec. IV A]. 

The construction described in this section should, in 
its simplicity, mainly be taken as a proof-of-concept il- 
lustration of how to meet the fundamental 'same error 
requirement' of Eq. (24). In practice, numerical opti- 
mization techniques may prove either helpful or indis- 
pensable, especially in the presence of additional control- 
input constraints. For example, if negative control inputs 
are not available to implement the required balance pair 
{IqiQ*) a-s (Q'Q,Qi/2), we can invoke a more general 
construction given in [,:)7]. The latter only assumes ac- 
cess to stretchable control profiles to implement {Iq,Q*), 
with a resulting sequence that is longer by a factor of 3/2 
with respect to {Q'Q, Qi/2)- The important point, how- 
ever, is that any DCG construction necessitates a recipe 
for the control inputs that will logically produce a non- 
trivial relationship among the corresponding EPGs, in 
order to evade the No-Go theorem of Sec. IIIC. 



C. Permissible controls for general linear 
decoherence and pure dephasing 

We next specialize the above general constructions 
to two physically relevant error models on n (driftless) 
qubits: (i) the generic linear decoherence model described 

by Eq. (5), where Q.^ = f^i^^ includes arbitrary single- 
qubit error operators, Eq. (6); (ii) the pure dephasing 
(or single-axis decoherence) model, which is formally ob- 
tained from the above generic case by assuming that cou- 
pling to B occurs along a known direction, say Z, in 



which case Ba'^ = for a = 1,2. We will refer to the 
resulting subspace of system-bath error operators. 



lie = span{Z(*) ® fiW} = VL 



z 

e 1 



(29) 



as the '(pure) phase errors'. The dephasing model is a 
common approximation often used in situations where a 
quantization axis defines the internal states of the indi- 
vidual qubits (such as in standard NMR QIP settings 
[!)]), and allows for significantly less resource- intensive 
FTQEC schemes to be devised in comparison to general 
linear decoherence [5 5]. This is not only due to the fact 
that the errors involved in dephasing are inherently re- 
stricted, but also to the fact that the algebraic structures 
associated with their growth are much simpler [51)]. 



1. One- and two-body controllable Hamiltomans 

Within the network model of QC, the majority of pro- 
posals invoke the use of control Hamiltonians that act 
non-trivially on single or pairs of qubits. Let us consider 
the set of primitive gates which are generated by switch- 
ing on and off control input in a control Hamiltonian 
parametrized as follows: 



H. 



(t) 



H. 



(pail-; 



Ctrl 



(30) 



(single 
Ctrl 



E 



E 

i=l a=X,Y 



(0 

a 1 



H, 



(pair) 



Ctrl 



i^j a=l 



where we generally allow for homogeneous two-body in- 
teractions which include natural entangling Hamiltoni- 
ans such as the Ising (ZZ) and Heisenberg (S • S = 
XX -\- YY -\- ZZ) interactions. We will focus on prim- 
itive gates during which the control propagator can be 
expanded as 



C/ctrl(t,0) 



exp 



n 3 



(31) 



i—l a—1 
i<j a—1 

where the functions Oa\t), 9a'^\t) are suitable control 
integrals and the set of qubits involved in two-qubit con- 
trol has no intersection with the set of qubits involved in 
single-qubit control. That is, we require that 

{z|0«(<)^O}n{j|e-''ni)^O} = 0, (32) 

so that A^i^{t) and A^2}{t) commute. For example, ap- 
plying a ZZ Hamiltonian between qubits 1 and 2, while 
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Figure 3: (color online) Left: Cayley graphs for Z2 (top) and 
Z2®Z2 (bottom), relevant to pure dephasing and arbitrary lin- 
ear decoherence, respectively. Representative Eulerian cycles 
and the corresponding EDD sequences are depicted. Right: 
Modified Cayley graphs for DCG constructions based on Z2 
(top) and Z2 ® Z2 (bottom), respectively. Relevant Eulerian 
paths and the corresponding DCG sequences are depicted. 



qubits 3, 4, and 5 are affected by single-qubit Hamiltoni- 
ans is an example where the above constraint is naturally 
obeyed. Note that such a 'divided control' is implicit in 
the circuit model of QC where gates are applied in paral- 
lel only when they commute. We emphasize that unless 
certain commutativity conditions are satisfied, in general 
Oa''\t) ^ ha''\s)ds, and a hnear relationship among 
the control inputs does not translate into a linear rela- 
tionship among the control integrals. 



2. Error-corrected gate constructions 



A good (minimal) DD group for the error set fll^^ 
corresponding to arbitrary linear decoherence is given by 

^Gi^jy^ 1.2(E) 1^2, 
I 



represented in Hs = (C^)®" through the n-fold tensor 
power representation {st^^^} = {Is,X'^^"\y'^^"\ Z'^^"^ 
[12, 35]. ^LD has two generators, which under the above 
representation we may take to be {Fj} = {X'-^^^\Y^'^^^^. 
Observe that the subspace of operators decoupled by Gld 
is not limited to as it includes, for example, inho- 

(i) (i) 

mogcneous bilinear couplings of the form Sa Sp \ with 



i ^ j and a ^ f3. Thus, flc,Li) ^ ^i^^ , which allows a 
wider class of primitive gates to be corrected using C/ld 
than the NOOP gate with single-qubit errors as originally 
analyzed in [■>'>]. 

The DD group for pure phase errors is simpler: 



Z2, 



represented via {Is, X'^^''^}, and with the single generator 
{F} = {X^*^'')}. Notice again that this representation 
also decouples errors beyond the simple dephasing terms, 
for instance inhomogeneous error operators of the form 

with i ^ j, therefore again we have Oc.z 3 rie ■ 
A pictorial view of the Cayley graphs relevant to both 
error models under examination is given in the left panels 
of Fig. 3, along with representative Eulerian cycles. 

A possible EDD sequence for correcting NOOP under 
the general decoherence model is given by: 



EDD'^ 



jj-(all)y(all)j|^(aU)y(all)y(aU)j^(all)y(all)j|^(an) 

(33) 

where operations are understood to be applied from 
right to left, following the underlying Eulerian cy- 
cle. This implies a total number of md = 
4(group elements) (2 generators) = 8 gates. Similarly, 
EDD under pure dephasing may be implemented by the 
sequence: 



EDD^ ^ (all) ^ (all) 



(34) 



which consists of just md = 2 gates, and may be viewed 
as implementing a form of 'continuous-time' spin-echo 
[57, 58]. The relevant DCG constructions are based on 
modifications of the above EDD gate sequences. Let Q 
denote a primitive gate. We define Iq = QQ' and Q1/2 
using the constructions of Sec. IV B. For general linear 
decoherence, a DCG for Q is then given by: 



(all) 



(35) 



r 



Similarly, for pure dephasing, we simply have: 
DCG^ ^ Qi/2X^^")/qA:(''"'. 



(36) 



The sequences given in Eqs. (35)-(36) represent Eulerian 
paths in the modified Cayley graphs corresponding the 



DD groups CJld and Qz, obtained as described in Sec. 
IV A; see also right panels in Fig. 3. If all primitive 
gates have the same duration r, since each added 'arm' 
in the modified graph has duration 2r, the resulting time 
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overheads may be summarized as follows: 





Linear decoherence 


Pure dephasing 


EDD 
DCG 


8 

16 


2 
fi 



3. Error per gate structure 

In order for the DQEC constructions just provided to 
be valid, it is necessary to explicitly show that the errors 
associated with the primitive gates generated from the 
(divided) one- and two-body Hamiltonians of Sec. IV C 1 
can indeed be decoupled by the appropriate DD group. 

Consider linear decoherence and ^Jld first. Let the 
control propagator for a primitive gate Q be given by 
Eq. (31). The toggling frame error Hamiltonian He{t) 
[Eq. (12)] is given by 

where we used the fact that [A^'^'^ (t) , A^'^'^ (t)] = 0. No- 
tice next that exp [i^'f^^(t)] iJg exp [— i^^^^(t)] can be 
expanded as a sum of single-qubit terms up to the first 
order Magnus, thus 

Up to the first order we may then write: 

^G.Mit)] = e^^'''(*)ng,jF,(t)]e-^^*''(*) = mods, 

where we used the fact that [6*"*'""*^*^ Sa^^^] = 0. Subse- 
quently, the first-order error $q is decoupled by Old- 

ne,o(*Q ) = ng,„ H,{s)ds>j = mods, 

which establishes the desired error cancellation. 

For the pure dephasing model, the set of allowed con- 
trol Hamiltonians in DCGs is slightly more restricted. 
For example, consider a DCG construction based on Eq. 
(36), where Gz is represented by {/, X^*^''^}. In this 
case, the single-qubit control Hamiltonians used in Q 
cannot include y*-'' terms as they would result in errors 
that are not decoupled by this particular representation. 
However, one can alternatively use the representation 
{/, oiQz, if the control Hamiltonians used in Q are 

given by 1"'*' . The choice of the particular representation 
of Q used in constructing DCGs may thus affect the set 
of permissible Hamiltonians. While this is not surprising 
in view of existing results on control of decoupled evolu- 
tions [-59], analyzing in full generality the compatibility 
of primitive gate sets with DQEC strategies for a given 
error model of interest warrants a separate investigation. 



V. CASE STUDY: DYNAMICAL CORRECTION 
OF SPIN-BATH DECOHERENCE 

In this section, we specialize the system-independent 
analytic DCG constructions described in Sec. IV to the 
paradigmatic setting of spin-bath decoherence. While 
quantitative modeling of a specific device is not our pur- 
pose here, our analysis is inspired by spin-based QIP ar- 
chitectures - in particular, electron spin qubits in semi- 
conductor quantum dots, as considered for instance in 
Rcfs. GO, 61, 62. In which physical parameter regime(s) is 
the improvement predicted for DQEC methods actually 
to be seen? What distinctive physical features and per- 
formance trade-offs are associated with DCG construc- 
tions? Beside complementing and expanding our previ- 
ous numerical results [23], these are two questions that 
the present investigation further addresses by example. 

A. Model system and primitive gates 

Let the system S consist of n individually addressable 
spin- 1/2 degrees of freedom, and let the environment B 
likewise consist of ti^ spin-1/2 particles. If we denote by 
I^*^ the spin vector operator associated with the i-th bath 
particle, the internal Hamiltonian Hint = Hb + Hsb we 
consider reads: 

Hb = E r,,(I«.lW-3/«4^'), 

l<i<j<nB 

Hsb = EgA^SW-lW. 

1=1 k=i 

Physically, Hb describes dipolar interactions between the 
bath spins, whereas Hsb describes a (hyperfine) contact 
interaction of each qubit with the bath. For concrete- 
ness, we further assume in what follows that the coupling 
strength Tij between bath spins i and j is arbitrarily 
chosen from a uniform random distribution in [— F, F], 

and similarly that the hyperfine coupling strengths A^'-' 
for each qubit are independently sampled uniformly at 
random from the [—A, A] interval. 

In line with Eq. (30), control over individual qubits 
and pairs of qubits is introduced through a gating Hamil- 
tonian of the form 

n n 
i=l a=X,Y = l 

where S*^'-* • S'^-'^ is the Heisenberg exchange interaction 
[(i3]. Thus, each primitive gate is realized by switching 

the control parameters h^^ (*) ' ^'y (0 ' ^'^'^ ^w''^ (0 ' under 
the 'divided control' assumption of Eq. (32). In particu- 
lar, we can simply assume that at each time t the qubits 
controlled through y(^), and S^^^ • S^-?) are distinct, 
and use as a starting point for DQEC the universal set of 
primitive gates given by ideal gates {exp(— i^C)}, where 
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C can be X^^\ y(^), and S(*) • S^^^. We restrict avail- 
able control inputs to the simplest choice of rectangular 
profiles: thus, a primitive gate is simply achieved by turn- 
ing on the corresponding control input for a duration r 
at a fixed strength h, allowing any control Hamiltonian 
to be written as a piece-wise constant function of time. 
The choice of the rectangular profiles allows for numer- 
ically exact simulations [(i4], while avoiding inessential 
complications in describing relevant control constraints 
and imperfections. In particular, the bounded-strength 
control constraint (cl) translates into demanding that 
maxt|ft,Q(t)| < /imax, while the constraint (c2) implies 
a shortest gating time r > T^mi [Sec. II B]. As an il- 
lustrative example of control errors, we may define over- 
rotation errors by replacing the nominal rotation angle 
26 in each (piece-wise constant) control segment with 
29{1 + e). If e fluctuates non-deterministically over dif- 
ferent instances of a particular gate, we have a random 
over-rotation; a constant (realization-independent) e de- 
fines a systematic over-rotation. 

Within the above setting, the internal Hamiltonian is 
the error Hamiltonian, H^. = ili„t and, as remarked, there 
is no drift, Hs = 0. Since Hsb belongs to f^l^^, as de- 
scribed in Sec. Ill the relevant DD group ^Jld discussed 
in IV C 2. Note that the DD problem for a single electron- 
spin qubit in contact with a nuclear-spin bath has been 
extensively analyzed using BB pulses in [G5, G6, 67], the 
limit Hs ~ corresponding to the low-bias regime where 
transverse and longitudinal relaxation fully compete. In 
our case, the two available group generators are realized 
through collective spin-flip gates, and the DCG construc- 
tions of Sec. IV A may be used. Thus, a primitive gate 
Q = exp(— i^C) of duration t is converted to a DCG gate 
through the following control Hamiltonian: 



B. Numerical results 

Perhaps the most important and obvious question, 
if DCGs are to be incorporated in a QC architecture, 
is whether they can usefully improve EPGs. While in 
theory DCG constructions are provably effective in the 
asymptotic limit of small errors, how stringent in practice 
this limit might be, can sensitively depend upon the spec- 
ifications and operating constraints of the device tech- 
nology at hand. For this reason, quantitative predictions 
on the effectiveness of DCGs are thus best formulated 
and addressed having a specific experimental platform in 
mind. Numerical simulations in toy models as we exam- 
ine can yet be instrumental in developing intuition on 
how to map out the actual regime of improvement as dif- 
ferent parameters are varied - in preparation for realistic 
scenarios where both Hsb and Hb arc typically fixed by 
either physical or fabrication constraints. 

In the numerical analysis that follows, we explore a 
wide range of open-system as well as control parameters, 
specifically A, F, and the base gate interval t, which for 
the current purpose can be also thought as identifying 
the shortest accessible modulation timescale, r = Tmin. 
In order to allow for arbitrary precision matrix calcula- 
tions (see also the Appendix), we focus on a relatively 
small open system, letting n = 2, ub = 6 henceforth. In 
all simulations, the qubit register is initialized in the fixed 
reference state |V'in) = (|00) -I- |01)/-\/2), whereas the en- 
vironment is initially maximally mixed, = Ib /2"^ . 
The evolution of the combined density matrix for the 
system and the environment is then calculated by mul- 
tiplying exact matrix exponentiation for piccewise con- 
stat Hamiltonians describing the segments of the evolu- 
tion. As representative target quantum gates, we con- 
sider single- and two-qubit rotations by 7r/4, that is: 



where explicitly we have (cf. also the DCG circuit given 
in Fig. 1 of [2:-i]): 



R 



H, 



+ec, 
-ec. 



-X(")), i = 1,7, 12,14, 
, i = 4,10, 11,13, 
i = 2,5,8, 
i = 3,6,9, 
i = 15,16. 



Notice that besides the original primitive gate Q used 
in segments i = 2, 5, 8, in segments i = 3, 6, 9 we imple- 
ment the gate (needed for the special NOOP gate Iq) 
through negative control inputs, while Q1/2 in segments 
i = 15,16 is obtained by using half the control power 
than for Q. 



(1,2) 
/4 



exp[-i(7r/8)X(i)](8)/(2), 
exp[-i (7r/8)S(i) • S^^)] = ^/i^^Xp, 



(1) 

■7r/4 



where the last equality makes explicit contact with the 
(universal) square-root-of-swap gate [2, 60]. Both the 
(first-order) improvement ratio introduced in Eqs. (2)- 
(3), r^^l = r, or directly the corresponding (in)fidelity 
contributions for the given initial state, will be used as a 
metric for quantifying DCG performance. 

We first focus on the case of perfect control (e = 0). 
Representative results on the dependence of the improve- 
ment ratio r upon the (minimum) gating interval r are 

depicted in Fig. 4, where the above W^y^ gate is con- 
sidered for a fixed realization of the dipolar-bath Hamil- 
tonian Hb [()N]. Weak- to strong- coupling dccohcr- 
ence regimes are then probed by varying the ratio A/T. 
All curves demonstrate the expected improvement in the 
limit of sufficiently small r. Different features are worth 
highlighting: first, as r 0, the data are consistent 
with the quadratic scaling with 1/t expected from Eq. 
(3) (for A = \, for instance, a fit of log(r) vs. log(T) 
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Figure 4: (color online) Improvement ratio as a function of 
(minimum) gating interval r, for the two-qubit target gate 
W^yf'' . The 'no improvement' (horizontal) line r = 1 is in- 
cluded for reference in this and all other figures that use r as 
a metric for performance. The internal-bath Hamiltonian is 
kept fixed at F = 1, whereas A/F takes the value lO'^ [brown 
diamonds], 10^ [green squares] and 1 [blue circles]. Inset: 
Infidelity error for the corresponding DCG implementation, 
1 — fwucG 1 ^ ^ function of r, for the same parameters used 
in the main panel. 



yields 1.9991 ±0.0002). Emergence of such a scaling may 
be taken as a good indication that asymptotic conditions 
arc well obeyed by the implementation. Second, the min- 
imum gate interval t = t*(A;T), at which the improve- 
ment region for DQEC (r > 1) is entered as t is de- 
creased, depends on the system-bath coupling strength - 
a larger A (stronger decoherence) requires shorter values 
of T* to be accessed in order for DQEC to be effective, as 
intuitively expected. Note that, qualitatively, a depen- 
dence of r upon both Hsb and Hb is expected on the 
basis of the uncorrected error associated with a DCG, 
Eq. (26). Only in the limit where F ^ (a so-called 
'non-dynamical' bath [G5, 67]) does the first-order im- 
provement regime simply relate to the coupling strength 
II -fossil I which is fixed through A. 

Exploring the interplay between Hsb and Hb as they 
are simultaneously varied for fixed (finite) values of the 
gating interval r reveals additional interesting structure, 
as shown in Fig. 5. In particular, the loss of fidelity 
of a primitive single-qubit gate (top panel) is contrasted 
to the loss of fidelity of the corresponding DCG version 
(bottom panel) over a range of intra-bath and system- 
bath couplings, F and A, respectively. Neither the un- 
corrected nor the DCG gate implementations depend in 
a similar manner upon the two parameters, the asym- 
metry in behavior being 'amplified' by the DCG. As a 
first interesting feature, only a mild dependence upon F 
is observed for primitive gates implementations, the cor- 
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Figure 5: (color online) Iso-fidelity contour plots for the im- 
plementation of the single-qubit target gate R^^j^ for fixed 
(minimum) gating time r = 1. Top: Primitive uncorrected 
implementation. Bottom: DCG implementation. In both 
cases, the integers labeling each curve give (the logarithm of) 
the corresponding infidelity error, with darker colors corre- 
sponding to higher fidelity values. 



responding fidelity increasing for very strong intra-bath 
dynamics. This effect can be partly understood if, as 
mentioned in Sec. IV A, Hb is included in the definition 
of the toggling frame: A very strong Hb will then in- 
duce fast oscillations in He{t), resulting in smaller EPGs 
[cf. Eq. (12)]. Secondly, and more interestingly, a non- 
monotonic behavior emerges for DCG implementations 
as the intra-bath coupling F is varied. 

The dependence of fidelities for DCG and primitive 
gates on Hb through F is further highlighted in Fig. 
6, where two representative values of t are examined. 
The resonance-like behavior in the DCG patterns is evi- 
dent, even when no significant improvement is expected. 
DQEC strategies such as DD or DCGs thus may probe 
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Figure 6: (color online) Infidelity, 1 — /, as a function of 
intra-bath/system-bath coupling scales for implementing the 
single-qubit gate R^^j^, at fixed A = 1. Primitive vs DCG 
implementations are contrasted for two different gating times 
T = 0.1 (top) and r = 0.01 (bottom). 
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Figure 7: (color online) Effect of systematic over-rotation er- 
rors in the DCG improvement ratio for Q — R'^^^ as a function 
of gating interval r. In both cases, the control error strength 
is parametrized by e € [0, 1], and we have set yl = 1, F = 0. 



the strength of the internal bath dynamics Hb, making 
them potentially attractive as a diagnostic tool for com- 
plex non-Markovian open quantum systems. 

Besides the error effects associated with Hsb and Hb, 
it is instructive to explore the impact of classical control 
errors in our setting, for example focusing on overrotation 
errors in the single-qubit gate R'^^J^- As further discussed 
in Sec. VIA below, two different models for control im- 
perfections may be physically relevant: (i) The nominal 
control Hamiltonian may be modified through the addi- 
tion of a 'deviation' Hamiltonian with a fixed strength, 
-^gatclO = -ffgate(i) + sH^cv, whcrc H^cv IS a Hermitian 
operator with units of energy and e a small number (re- 
call that Hgate{t) = Hctri{t) as long as Hs — 0). This 
results in a scaled systematic EPG, that can be tolerated 
(to first order) DCG constructions, similar to the errors 
due to the bath, (ii) The nominal control Hamiltonian 
may be systematically scaled by a factor: 

^^gatcW = (l+^)^gatcW- (38) 

This results in a fixed systematic EPG, in the sense that 
even in the absence of the bath, a control segment of 
duration r intended to produce exp{—i9C) will instead 
implement exp[—i9{l+e)C], regardless of r. In the simu- 
lations, model (ii) is implemented, given that there is no 
a priori indication that such errors should be tolerated 
by DCGs, unlike errors of type (i). 

Illustrative results for the DCG improvement ratio r 
obtained for a non-dynamical bath (T — 0) at fixed 
A = I are summarized in Fig. 7, as a function of the 
gating interval and for different values of the control er- 
ror strength. Overall, although the observed performance 
remains fairly stable for error strengths up to a few per- 
cents, no intrinsic fault-tolerance is observed for arbitrary 



strengths and/or rotation angles. In particular, no im- 
provement is achieved if uncorrected control errors dom- 
inate over the decoherence errors (due to the spin-bath 
in this case) that DGCs remove, as anticipated in [23]. 
For small values of e (roughly below 1%), there is a set of 
values of r for which the decoherence error is indeed the 
leading source of gate error, and for this region r consis- 
tently exceeds 1. For the remaining values of r, and/or 
too large values of e, no improvement is observed: either 
T is too large for the asymptotic DCG regime to set in 
(cf. the interval [r*, oo] in Fig. 4, as discussed earlier), or 
T is too small and the decoherence errors become negli- 
gible compared to the uncorrected control errors. In the 
latter case, no benefit arises from further reducing t, and 
performance saturates at an e-dependent limiting value. 
Qualitatively similar results may be obtained for random 
control errors. 



VI. EXTENSION TO FAULTY CONTROLS AND 
SYSTEMS WITH DRIFT 

The explicit constructions of Sec. IV for DCGs rely 
on a straightforward relationship between the intended 
gate Q and the EPG associated with it. This is en- 
abled by assumptions (al) and (a2) in Sec. HA: we 
have focused on Hsb (causing decoherence) as the sole 
source of (decoherence) errors, and the gating mecha- 
nism neither employed nor was corrupted by drift terms 
in Hs- In reality, as mentioned, classical errors due to 
control imperfections and additional quantum errors due 
to non-trivial unitary evolution are likely to be present 
to a lesser or greater extent. As also mentioned, robust 
control methods have been devised for separately tack- 
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ling these errors, whose impact may be deleterious even 
if isolation from the environment is good hence decoher- 
encc is not a primary concern. In a way, composite pulse 
techniques for systematic control errors [25, 2G, 27, G9] 
may themselves be considered as a classical dynamical er- 
ror correction strategy: they are pcrturbativc in nature; 
they assume no quantitative knowledge of the value of 
errors and arc constructed from building blocks of primi- 
tive pulses. Likewise, strongly modulating and gradient- 
ascent pulse engineering techniques in NMR [29, 30, 31] 
share with DQEC the basic idea of suppressing unwanted 
(but known) coupled-spin dynamics through coherent 
time-averaging, numerical search being explicitly invoked 
to synthesize an optimal modulation. 

The above considerations prompt the need for explor- 
ing DQEC under more relaxed control assumptions than 
(al) and (a2). In particular, a natural question is 
whether it is possible to produce hybrid dynamical error 
correction strategies to simultaneously counter different 
types of classical and quantum errors, over which a differ- 
ent degree of knowdlege may be available to begin with. 
While this question is too broad in scope to be answered 
in full here, and might be best treated in conjunction 
with numerical optimization approaches, wc nonetheless 
illustrate the flexibility of analytic DQEC methods by 
providing three specific examples of DCG DCG construc- 
tions beyond (al) and (a2). 



Towards hybrid composite-pulse/ 
DCG constructions 



A large class of systematic control errors may be natu- 
rally modeled at the Hamiltonian level. That is, recalling 
Eq. (8), we may assume that the nominal gating Hamil- 
tonian ffgate(i) during the interval [ti,t2] is modified as 



-f^gato(i) — -f^gato(i) + Hs,g + sHdcv{t), 



(39) 



where i?dev(i) is (typically gate-dependent) Hamiltonian 
whose functional form is known, and e is (possibly un- 
known) real number such that 



C^gate(^2j ^l) ^ 7 + 



exp 



I I Hcti-iis)ds 



< 1. 



The systematic nature of the error is implicit in the ex- 
istence of a functional dependence between i?gate(0 and 
Hiiev{t) as a mapping between time-dependent functions. 
For physically relevant error models, this mapping is of- 
ten specified in terms of simple relationships between the 
control and the error Hamiltonian - e.g., for a single 
qubit, application of an X gating Hamiltonian may in- 
duce a 'parallel' (in-axis) deviation Hamiltonian also pro- 
portional to X, or a 'perpendicular' (off-axis) deviation 
Hamiltonian proportional to Y,Z (see also [49, G7]). 

More generally, a relatively simple description of sys- 
tematic control errors may be established by directly in- 
corporating and examining their effects into the error 



action operator. Thus for each target gate Q ideally 
implemented via ifgate(i) but actually implemented via 
'^gatc(^) in the presence of the bath, to first order in the 
gate duration we can expand the associated total error 
action $ as: 



$ = $SB + ^'ctrl, 



(40) 



where for simplicity we have assumed that no contribu- 
tion arises from Hs (or that the latter commutes with 
HsB and $ctri accounts for systematic control errors. 
The assumptions (al)-(a2) allowed us to let $ctri = 0, 
and provided indirect control over ^sb through the gat- 
ing interval r: For primitive gates in the absence of non- 
commutative drift terms, ^sb is linearly proportional to 
T. Such a dependence of ^sb on t could then be explic- 
itly utilized in the constructions of Sec. IV B for match- 
ing the errors in QQ' and Qi/2. 

Fixed systematic errors [cf. Eq. (38)] refer to <f>ctri 
terms in Eq. (40) that do not scale with t, as opposed 
to scaled systematic errors that scale linearly with t. 
Hamiltonian-control errors as described by Eq. (39) al- 
ways lead to scaled systematic errors. Since the latter 
can be formally absorbed into $sb, in the absence of 
additional (drift) errors, the constructions of Sec. IV B 
may still be used to effectively reduce the effect of sys- 
tematic control errors in the same way the effects of the 
bath are corrected. That is, arbitrary DCGs are intrinsi- 
cally robust against scaled systematic errors. The FDD 
construction for NOOP on the other hand (as a special 
DCG), is tolerant of all, scaled and not, systematic con- 
trol errors (as long as they are small), for they as can be 
absorbed into generator gate errors p. [35] . The numer- 
ical results of Sec. VB confirm instead the expectation 
that generic DCG gates need not be robust under fixed 
systematic errors. 

It is still plausible that one may use a combination of 
composite pulses and DCGs to actively stabilize quan- 
tum gates against both dccohcrence and systematic er- 
rors. The starting point is a composite pulse construction 
that in the absence of the environment can generate a de- 
sired gate Q robustly with respect to systematic errors 
within a specific model. Let us assume that it is possi- 
ble to modify the control inputs for this construction to 
obtain the gates Q1/2 and Q defined in Sec. IV, such 
that Qi/2 and Q' constructed in this way are also com- 
posite pulse robust against the same control errors. Then 
clearly the DCG construction that employs Q, Q' , Q1/2, 
along with the generator gates Fi , will correct (to first or- 
der) the systematic errors and the errors due to the bath. 
This, however, necessarily involves longer and more com- 
plex sequences of primitive control Hamiltonians. 



B. Qubits in rotating frames 

An important class of systems with drift are qubits 
which are defined in terms of two (non-degenerate) in- 
ternal energy levels and are addressed using resonance 
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techniques. This is, for instance, the typical setting in 
both hquid- and sohd- state NMR QIP devices [9] . In this 
case, the internal open-system Hamiltonian Hint includes 
'chemical-shift' terms of the form Hs = J2i ^i^i: which 
are employed for single-qubit manipulations by moving to 
a frame which rotates with the applied carrier frequency 
flc and by tuning on-resonance with the target spin(s). 
The rotating frame transformation has the effect of in- 
troducing an explicit time-dependence in system (hence 
system-bath) operators, according to 

O d{t, to) = cxp{~itnc Zi) O cxp{imc ^ Z,), 

i i 

where we have assumed that to = 0. Thus, even ne- 
glecting contributions to Hs from off-resonance qubits, 
both the interaction with the environment and therefore 
the error Hamiltonian arc transformed, in general, into 
carrier-modulated Hamiltonians HssitjO) and He{t,0), 
respectively. This causes most of the constructions in 
the current paper (including EDD) to be not immedi- 
ately applicable to compensate the rotating-frame error 
associated with a gate, because this error will generally 
depend on the time at which the gate is applied. Still, 
thanks to the known periodicity of the time dependence 
of the error actions, it is possible to synchronize the du- 
ration of the applied primitive gates with the rotating 
frame period 27r/fi, in such a way that the relationships 
between the errors assumed in EDD and DCG construc- 
tions are stroboscopically preserved. 

Furthermore, in practice the frequency Qc may be 
much higher than typical inverse gate durations. This 
results in rapidly oscillating terms in the rotating francs 
that are effectively canceled [when Eq. (11) is used with 
He{t)] . If He consists of arbitrary single-qubit error terms 
[Eq. (5)], then in the limit of high carrier frequencies, 
the effective error model for the open-system dynam- 
ics in the rotating frame approaches the pure dephasing 
error model, since any terms proportional to X'*^ and 
y^*-* in the error action will be averaged to zero, whereas 
any terms proportional to Z^*-* will be time- independent. 
This has two implications: On the one hand, single-qubit 
drift terms in Hs are of no concern for DQEC in the im- 
portant case where the physical error model is dephasing 
to begin with; on the other hand, if arbitrary decoher- 
ence is the appropriate model in the physical frame, DCG 
constructions for pure dephasing are directly applicable 
in the rotating frame in the limit of large f2c- 



C. Always-on qubit-qubit interactions 

Aside from single-qubit drift terms as discussed above, 
another common type of drift contributions may arise 
from always-on two-qubit couplings in Hs ■ In NMR QIP, 
for instance, the latter can either be 'J-coupling' Ising- 
like Hamiltonians for weakly interacting spin systems, 
or Hcisenbcrg-like Hamiltonians in the strong-coupling 



regimes. Within our present approach, the main com- 
plication that a non-switchable component Hs^g in the 
gating Hamiltonian for Q introduces towards construct- 
ing a DCG version is in establishing the existence of 
the required Q1/2 and Q' gates, since Hs^g effectively 
imposes a limit on the range of achievable control in- 
puts. Specifically, the constructions of Sec. IV depend 
on generating propagators Ug^^^{t, 0) and U'^^^^{t, 0), and 
it is not immediately clear whether universal control over 
the system S suffices for realizing such propagators (as 
a one-parameter family of time-dependent operators as 
opposed to just end-point unitaries), even in an approxi- 
mate sense. Nonetheless, we see no fundamental obstacle 
to treating such cases through strong modulation of the 
system's dynamics and numerical optimization. To fur- 
ther support our optimism, we show here how analytical 
DCG constructions remain feasible upon careful consid- 
eration of the relevant constraints on a case- by-case basis. 

Consider a system of n Hcisenberg-coupled qubits in a 
nearest-neighbor configuration: 

n-l 

i/s = A^S«-S(*+i\ (41) 

1=1 

where A G M is a fixed coupling strength. Given Hs 
in Eq. (41), universal control over S can be gained by 
assuming access to a switchable control Hamiltonians of 
the following form: 

H.try{t) = 5^[4>(i)A« + (42) 

i 

where as usual the control inputs h'^a {t) arc subject to the 
power and bandwidth constraints (cl)-(c2) . We assume 
the bath to interact with S according to the arbitrary 
linear decoherence model: Hsb G ^e^^ ■ Our goal is to 
provide a (first-order) DCG construction for the following 
gates: exp[— i6'S'i*'''], where a — X^Y, and cxp[— i6'S('^^ • 
g(fe+i)j^ such that the resulting EPG scales with t'^^^- 
Our methodology here will be somewhat different from 
Sec. IV, in the sense that primitive gates are not directly 
used as the basis of DCG constructions. 

We first make the following general observation: Sup- 
pose that the gating propagator and the error Hamil- 
tonian for a control block in the interval [0, T] can be 
respectively written as [/gato(0,t) ~ Ug^i{t)Ug^2{t) and 
He = He,i + He.2, in such a way that 

[UgAt),UgMt')]^0, (43) 

[Ug,i{t),He,2] = [Ug^2{t),He,l] = 0. (44) 

Using the above equations and the tools of Sec. II C, we 
can write the error associated with this block up to the 
first order as <I>[^1 = -I- $2^', where 

= / Ug.iity He,iUg,iit)dt, (45) 
Jo 

4'' = / Ug,2{t)^He,2Ug,2{t)dt. (46) 

Jo 
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As it turns out, generating a two-qubit gate on a target 
qubit pairs, exp[— i^S^*^^ • S^*^"^^)], is the most non-trivial 
step. In this case, we clearly have Hs^g — S^'^^ • S''^+^' 
and Hs — Hs^g = Hs,e- Basically, our strategy is to 
first devise a scheme for dynamically correcting all error 
terms affecting any qubit other than k and fc + 1, and 
next to superpose a scheme that remove any additional 
error term while steering the target qubits according to 
Hs,g. For the first step, consider the DD group C/nn = 
Z2 (2) Z2 (81 Z2 (8) Z2 , designated to decouple the subspace 
of bilinear qubit terms 



with 



St 



where 



1, • • • , fc - 1, A: + 1, • • • , 71 - 1, Si'^ e B{Hb) 
and Sp'^^'^ may include the identity operator C/nn 
has 16 elements, and 4 generators which we choose to 
represent in Ti.s via 



where 



(47) 



;^odd 



ncU) ccv _ TT cU) 
'-'a T '-'a ^ _[_!_ '-'a i 



Ek = {j I j even; < j < n; j 7^ fc, fc + 1}, 
Ok = {j I jodd; < j < 7i; j^k,k+l}. 

Consider the EDD sequence based on {-Pi.i}, and let 
Hg^i{t) be the gating Hamiltonian associated with it over 
the interval [0, Ttwo] = [0, 16 x At], assuming as usual that 
each generator is implemented over an interval of length 
r. Note that Hgj{t) can be generated through appropri- 
ate combinations of the switchable inputs in Eq. (42). 

In the second step, we construct the entangling gate 
exp[— 64iTAS('^' • S^'"'"'"^^] by applying another control 
Hamiltonian Hg^2{t) in parallel with Hg^i{t). Hg^2{t) im- 
plements a sequence very similar to the EDD sequence 
for the group ^Jld = Z2 (8) Z2 designated to cancel errors 
in the subspace 



LD.k 



{5i*)®sW|i = fc,fc + l}. 



However, while the original EDD sequence would employ 
the generators X^^'^X^^+^'^ and ( or a similar 

combination), we now replace the two original bilinear 
generators with 

Fi,2 = xWx('=+i)exphl6iTASW -SC^+i)], 
F2,2 = yWy('=+i)cxp[-167rAS(^-) -8(^+1)]. 

The modified EDD sequence can be implemented through 
the available control Hamiltonians via a collective gating 
Hamiltonian of the form 

Hga{t) ^ hI''\ for t e [8{l - l)r, 8/t], 



^(fe) ^ I Tfr / = 1, 3, 6, 8, 
' |^[F('=)+y('^+i)]+Hs,«, 1 = 2,4,5,7. 

Notice that the system gating Hamiltonian, explicitly in- 
cluded in Hg^2{t), commutes with the single qubit terms 
acting on fc and fc -|- 1 in Eq. (48). Let 

n n 



1=1 

i=jLk 



i=l 
i^k,k+l 



k+1 



a i—k 



We may thus readily verify Eqs. (43)- (44). Also no- 
tice that the total gating Hamiltonian commutes with 
g(fe) . gC^+i). Thus, the overall action of the con- 
trol block defined via Hg^i{t) -\- Hg^2{t) is given by 
C/g,i(64T)t/g,2(64r) = exp[-64iTAS(''^) • and we 

can derive the error associated with the whole sequence 
(of length Ttwo = 64t) by invoking Eqs. (45)- (46): 

- rUgAtyH,sUgAt)dt = J2'^i\' 
Jo I 

4" - / C/9,2(i)ti^e,2f/g,2(t)rft-E4H> 

Jo , 



where and respectively denote the errors due to 
He,i and He,2, associated with interval [{I — 1)t,It]. By 
(EDD) construction 



[1] 



4 

E 



ne™($F, J. 



One can verify directly [similar to Sec. IV C 1] that 
using Hg^i{t) results in an error ^Fij associated with 

|2| 

each generator Fij that belongs to fl^. , and thus we 
have $ 



[1] 



mods. Similarly, one may also verify 
that the error $2^J belongs to fli,D,k- We may then use 
Eq. (15) and the fact that the commutator between 



ASC^) • SC^+i) and J^'r+l S^^^ ( 
to 0[r^ max( II II) A] to write 



Ena.B 



i?Q^ can be ignored up 



therefore (f>2^' = modi?. The total error associated with 
this sequence is thus modulo pure bath terms and terms 
that scale with r^, as desired. 



19 



Implementing a DCG version of a single-qubit ro- 
tation Q on a selected target qubit k is conceptually 
more straightforward. In this case, Hs = Hs\e thus 
He = HsB + Hs, and the DCG construction of Sec. IV 
can be reproduced with ease by letting the relevant error 
subspace cover all nearest-neighbor bi-linear interactions 
of the form 

where i = 1, • • • , n— 1 and Ba'^ e B{Hb)- The space fll'^^ 
can be decoupled by the group CJnn with a representation 
similar to Eq. (47), except that now the sets Ek and Ok 
are chosen as 

Ek = {j I j even; <j <n}, 
Ok = {j\j odd;0 < j < n}. 

Similarly, we may generate Q, Q', and Q1/2 by simply 
switching the appropriate single-qubit Hamiltonian on 
the k-th qubit. One may verify that the error associ- 
ated with each control segment belongs to fll . The 
final DCG construction will thus contain 64 segments for 
EDD and 32 segments for Q1/2 and QQ' , for a total du- 
ration of Tsingio = 96r. 

It is worth noting that, in a way. the larger number 
of control segments emerging from the above analysis as 
compared to driftless scenarios may be simply regarded 
as an analytically generated 'digitized' pulse shape, with 
a basic time-step determined by T-^i-a- Attention must be 
paid, however, to the fact that longer control sequences 
might potentially invalidate the (sufRcient) convergence 
assumptions implicit in DCG constructions. This high- 
lights the importance of finding an optimal DCG con- 
struction for a given control scenario. In this sense, the 
constructions of this section should be considered as a 
proof-of-concept as opposed to optimal results that might 
best be tailored to specific control scenarios and will likely 
involve a combination of analytical and numerically op- 
timized control designs. 

VII. CONCLUSION AND OUTLOOK 

We have provided the mathematical and control- 
theoretic framework for DCGs as introduced in Ref. 23, 
and further explored their range of applicability. Besides 
examining in full detail the analytical construction of 
decoherence-protected gates in the prototypical scenario 
of arbitrary single-qubit errors, a physically motivated 
setting involving spin-bath decoherence has been char- 
acterized through exact numerical simulation. This has 
both confirmed the expected performance of DCGs in a 
wide parameter regime and yielded insight into the inter- 
play of different open-system and control features. Addi- 
tional challenges arising in the synthesis of DGCs under 
realistic control assumptions involving drifts and control 
errors have been laid out, and proof-of-concept solutions 
constructed in specific representative situations. 



Our present analysis naturally points to several direc- 
tions for further investigation. First, our focus here has 
been on DCG constructions that cancel errors only up 
to the first order in Magnus expansion. In principle, it 
is possible to go beyond the first order DCG construc- 
tions and produce even more efficient DQEC protocols 
[37]. While this will unavoidably involve a higher level of 
sophistication in the resulting control modulation and/or 
the computational effort involved in numerical search 
processes, a clear understanding of relevant complexity 
properties is needed. Making contact with the concept 
of 'quantum control landspace' as recently proposed for 
the generation of unitary transformations in closed quan- 
tum systems [70] might prove especially insightful. 

Second, we only considered a limited set of permissible 
control Hamiltonians for DCG constructions. A general 
analysis of what we referred to as the 'gate permissibility 
problem' for a given open quantum systems is, as men- 
tioned, still lacking. This generalization appears even 
more interesting by observing that the set of permissible 
controls can be directly related to the specific representa- 
tion of the DD group used in the DCG construction. In 
fact, the apparent freedom in the choice of DD represen- 
tation independently motivates the existence of DQEC 
proposals for the same tasks as DCGs that do not rely 
upon DD explicitly. Numerical optimization is, in this 
sense, a compelling option to consider for investigating 
more general DQEC constructions. 

Lastly, different options exist for incorporating DQEC 
protocols into actual QIP architectures. While we have 
stressed the potential usefulness of DCGs directly at the 
physical level, DCG constructions could potentially be 
envisioned on top of a logical encoding, for instance to re- 
duce EPGs in silicon double-quantum-dot logical qubits 
as considered in [71]. Once an algorithm is specified, 
whether at the physical or logical level, one may further 
consider whether to produce the DCG version of each cir- 
cuit gate, and obtain an error-suppressed algorithm ac- 
cordingly; or whether instead to dynamically correct the 
whole algorithm viewed as a complex unitary transfor- 
mation. While a thorough understanding of such design 
issues is fundamental for real-world applications, optimal 
solutions will ultimately depend on various details and 
constraints in place and, as such, are most effectively ad- 
dressed in specific contexts. It is our hope that this work 
will prompt further interest and investigation in DQEC 
approaches from the QIP implementation community. 
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Appendix A: NUMERICAL METHODS 

The simulations described in Sec. VB were performed 
(in Mathematica) using direct exact exponentiation for 
obtaining the propagators describing the system-plus- 
bath dynamics in the presence of control. Arbitrary pre- 



cision matrices turned out to be necessary for the follow- 
ing reasons: (i) Small errors correspond to toggling-frame 
error propagators that are close to identity, thus it is es- 
sential that any numerical error be much smaller than the 
physical effect being probed, (ii) Each matrix multipli- 
cation effectively reduces the precision of the calculation 
by a fixed amount which depends on the dimension of the 
matrices used. We were limited to consider a relatively 
small system/bath and use matrices of very high preci- 
sion [60 digits or higher] throughout the computations. 
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